Exact matrix product solution for the boundary-driven Lindblad XXZ-chain 
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We demonstrate that the exact non-equilibrium steady state of the one-dimensional Heisenberg XXZ spin 
chain driven by boundary Lindblad operators can be constructed explicitly with a matrix product ansatz for the 
non-equilibrium density matrix where the matrices satisfy a quadratic algebra. This algebra turns out to be 
related to the quantum algebra U q [SU(2)]. Coherent state techniques are introduced for the exact solution of 
the isotropic Heisenberg chain with and without quantum boundary fields and Lindblad terms that correspond 
to two different completely polarized boundary states. We show that this boundary twist leads to non-vanishing 
stationary currents of all spin components. Our results suggest that the matrix product ansatz can be extended 
to more general quantum systems kept far from equilibrium by Lindblad boundary terms. 
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The non-equilibrium behaviour of open quantum systems 
has become accessible through recent advances in artificially 
assembled nanomagnets consisting of just a few atoms QJJ] or 
in the study of quasi one-dimensional spin chain materials like 
SrCuC>2 where many transport characteristics are measurable 
experimentally JUS]. In particular, it is desirable to under- 
stand the interplay between many-body bulk properties (e.g. 
magnon excitations or magnetization currents in quantum spin 
systems) and local pumping (applied to the boundary of a sys- 
tem) driving the system constantly out of equilibrium. A good 
starting point is provided by the anisotropic Heisenberg model 
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of coupled spins. The pure quantum version of this model 
is exactly solvable by the Bethe ansatz. Interestingly, within 
linear response theory, i.e., close to equilibrium, it was found 
that at finite temperature a diffusive contribution to the Drude 
weight appears 15|-[3], which is at variance with the long-held 
belief that integrability protects the ballistic nature of trans- 
port phenomena. Unfortunately the Bethe-ansatz fails in the 
more relevant context of open far-from-equilibrium systems 
where these questions can be addressed directly in terms of 
the Lindblad Master equation [8] 
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[H,p]+V L (p)+V R (p) 
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for the reduced density matrix p associated to the chain (here 
and below we set h = 1). The dissipative terms V L,R (p) = 
D L,R pD L,R\ _ i j p) d l,r d l,r1 I with the Lindblad opera- 
tors D L ' R acting locally at the open ends of the quantum chain 
(see below) describe the coupling to external reservoirs that 
drive a current through the system and thus keep the system 
in a permanent non-equilibrium steady state. Indeed, using 
dissipative dynamics for the preparation of quantum states is 
becoming a promising field of research l^ fioll . 



Significant progress has been achieved very recently in two 
remarkable papers by Prosen ITU [l2tl who observed that the 
exact stationary density matrix for the XXZ chain with one 
specific pair of Lindblad boundary terms can be constructed 
explicitly in matrix product operator form lfl3ll by a matrix 
product ansatz (MPA) somewhat reminiscent of the matrix 
product ansatz of Derrida et al. ITill for the stationary dis- 
tribution of purely classical stochastic dynamics. With an 
explicit representation of the matrix algebra Prosen was then 
able to compute analytically various physical quantities of in- 
terest. However, in contrast to Q, where the matrices satisfy 
a quadratic algebra, the matrices of [fTTl [l2tl satisfy a cubic 
algebra which arises from a peculiar local cancellation mech- 
anism involving three neighboring sites in the quantum chain. 
This feature is significant since, due to the lack of a general 
representation theory for cubic algebras, this approach does 
not lend itself easily to generalization to other open quantum 
systems with other cubic algebras or even small modifications 
of the original problem such as boundary fields or other Lind- 
blad terms for the XXZ chain which would require a different 
representation. Indeed, the wide applicability of the MPA of 
IU4I1 derives from the fact that many quadratic algebras (which 
include all Lie algebras through their commutation relations) 
have explicitly known representations which is crucial for the 
exact computation of physical observables lfl5tl . 

In this Letter we show that exact non-equilibrium steady 
states for open quantum systems can be obtained from a ma- 
trix product ansatz which yields a quadratic algebra. Specif- 
ically, we consider the Lindblad quantum XXZ chain and 
show that the associated matrix algebra is related to the bulk 
symmetry of the XXZ-chain, which is the quantum algebra 
U q [SU(2)] with A = (q + q~ 1 )/2. A coherent state rep- 
resentation makes it possible to consider Lindblad operators 
that correspond to two different completely polarized bound- 
ary states, viz., in the (y, z) plane on the left boundary 
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and in the (x, z) plane on the right boundary with 



(cos 8ro% - ia v N + $mQ R (j z N ). 
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These Lindblad generators lead to local dissipative terms 
whose stationary solutions, satisfying T> Li - R:> (p L ( R ') = 0, 
are respectively the pure states p L = 1/2(1 + er*) = 
ItuXttil and p R = 1/2(1 - cr*) = where 
) is the eigenstate associated to the eigenvalue +1 of 
er* = s'm0L<j y + cos#l<t z , and \i v ) is the eigenstate of 
a z = — sm9R(j x + cos 9r(t z with eigenvalue —1. For com- 
putational convenience we have chosen equal left and right 
amplitudes T in (O and By a judiciously chosen simi- 
larity transformation these amplitudes can be made different 

Moreover, we allow for quantum boundary fields acting on 
the directions of the local polarizations specified by (0 and 
(3J. Consequently we add to the Hamiltonian (Q~|i the con- 
tribution f L ■ a = f L o- z l for the left-end of the chain and 
f R ■ a = f R a z coming from the right-end boundary field. 
For convenience we choose J = 1/2 and e = 1 so that 

H = X^^Ti 1 h-k.k+i + 9i + 9n with the four-by-four ma- 
trix h = \{a x ® o x + (jy ® (jy + A(a z <g) a z - 1) for the 
nearest neighbour bulk interaction and the two-by-two matri- 
ces g L = f L a z and g R = f R a z for the boundary fields. The 
subscript indicates on which sites of the chain the quantum 
operators g and h act non-trivially. We write the stationary 
density matrix satisfying 



i[H,p]=V L (p)+V R (p) 



(5) 



in the standard form p = SS*/ Tr(SS^). 

Our starting point for solving d5j is a matrix product ansatz 
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which we augment by auxiliary matrices 5 such that the local 
divergence condition 



[h, f2(g)fi] =S®f2-f2®S 



is satisfied. In this construction 
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are two-by-two matrices whose entries A a ,E a are non- 
commuting matrices that act in a space A with inner product 
( • | •), and \ip) are vectors in A. In terms of Pauli ma- 
trices cr ± = \(o x ± i<J v ), cr 2 , and the two-dimensional unit 

matrix 1 one can conveniently write O = Aq\ + A ■ a with 
A = \{A t + A 2 ), A z = \{A X - A 2 ), A x = \(A+ + A.), 
A y = t;( A + — In our construction the local divergence 
condition leads to a set of 16 quadratic relations for the eight 
matrices A a , E a and the problem to be attacked is the con- 
struction of matrices which satisfy these relations. 

Remarkably, all 16 equations (0 can be solved in terms 
of only three independent matrices A±, Q with QQ^ 1 = 



Q 1 Q = 1 by choosing the auxiliary matrices E± 

(q - q-^bQ - cQ- 1 )^, E 2 = -(q- q^){bQ - 
setting 



A l =bQ + cQ~ l 1 A 2 = bQ + cQ~ 



= 0,E 1 = 
cQ- 1 )/2, 



(9) 



and requiring 



[A+,A-] = -(q - q-^ibbQ 2 - ccQ- 2 ) (10) 

(11) 



QA± = q ±l A ± Q 



The constants b, c, b, c are arbitrary. Choosing the 
parametrization 
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and defining 



then leads to 



A± =: iaS±, Q =: Xq * 
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These are the defining relations of U q [SU(2)], the q-deformed 
universal enveloping algebra of the Lie algebra SU (2), which 
is the non-abelian symmetry of the bulk Hamiltonian (Q} lfl7ll . 

After deriving a matrix algebra from the bulk interactions 
the second step is the explicit construction of such matrices 
and of the vectors ( V | and | W ) using the boundary inter- 
actions. For the present case we note that the representa- 
tion theory of U q [SU(2)} is well-understood and analogous 
to that of SU (2), except when q is a root of unity, where some 
special features arise lfl7ll . In particular, with the definition 
M — i -i — we nave the irreducible representation (irrep) 
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where p is an arbitrary complex parameter. This irrep is 
infinite-dimensional, except when 2p € N). The bra's and 
ket's form an orthogonal basis of A — C N with inner product 
( k | k' ) = 5k,k'- Relations (O, ([Pil l then provide a represen- 
tation of the matrices A a . 

In order to satisfy the boundary conditions involving the 
quantum boundary fields and the Lindblad dissipators we 
further define $ = [g , fl] and the tensor products = 
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The local divergence condition implies [ H , tt® N } - 
Sx + — Stv- The stationary Lindblad equation Q 
be split into two equations 

V L (S&) = + Hi)S t -i5 , ($f' + + 5|) 



(17) 
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for each boundary. Using the decomposition S = (<j>\ [I-Aq + 
cr z A 2 + A + cr+ + A_o~ )] ® ft®^- 1 1 ip ) for the first equation 
and S = ( If^" 1 ® [A 1 + A z a z + A+a+ + A^a~] \ i/> ) 
for the second equation and factoring out the term containing 
j}®jv-i yields ( WO se ( S f equations for the action of the ma- 
trices A a on the vectors ( <f> | and | ip ) respectively. In this let- 
ter we outline this programme for the isotropic chain A = 1. 
The construction for A ^ 1 is conceptually similar, but tech- 
nically more involved and will be presented in a more detailed 
paper Qj|. 

For taking the isotropic limit q — > 1 we choose the normal- 
ization factors a = A = 1 and set fi = v in (TT~2T >. ( fT3l l and 
arrive at 
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iS- v~ 1 S z 
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The quadratic relations for the quantum algebra turn into the 
usual commutation relations [S + , S-] = 2S Z , [ S z , S±] = 
±S± for SU(2). The irreducible representation ( fTTI i turns 
into an irrep of SU(2) by observing that [a;]i = x. Since the 
Lindblad dissipators do not generate terms proportional to the 
unit matrix, we cancel these terms that appear on the r.h.s. of 
(0 by setting v = i which leads to Aq = and ft = iS ■ a, 
where S = ( S+ ~^- ,i S+ ~ S - ,S Z ) and a = (a x ,cry,a z ), or 
in terms of the a the form ft = i(S z a z + S+a + + S-cr~). 

The key step in solving the boundary equations is the intro- 
duction of coherent states 
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Using the commutation relations of SU (2) one finds 

(<t>\S z = (<t>\(p-cf>S + ), {cf>\S-=ct>(0\(2p- ( j ) S + ) (22) 
and 

S z \iI>} = (p-tPS-)\tP), S + ^)=V(2p-^_)|V) 

(23) 

The left boundary equations can now be solved by noting 
that the Lindblad operator (0 can be obtained from a com- 
plete polarization along the z-axis by the unitary transforma- 
tion U = e l ^r a on site 1 of the chain which rotates the z-axis 
into a new direction u. After this transformation the leftmost 
matrix ft in the tensor product ft® N reads in the new basis 



with the new components 

S Z (6 L ) = S z cosO L +ismO L 

s+(o L ) = s+ + s - +cose L 

S-(0l) = S+ + S - -cos6 L 
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iS z sin6»j(25) 
iS z sin 9l- 



In order to solve the left boundary equation we need to impose 

(<t>\S_(9 L ) = 0, {ct>\S z {e L )=p{ct>\ (26) 

Using d22l these conditions are satisfied if the coherent state 
parameter <j> is chosen to be 



4> = itan(0£/2). 



(27) 



In order to prove this result we point out d22l . ( 1231 can be 
used to express vectors of the form ( <j> | (a+ bS z + cS+ + dS- ) 
that appear in the boundary equations just in terms of e.g. 
(0| (a' + d'S-), and similarly for the action on ket-vectors 
| ip ) . The choice (127] ) leads to 

S = {(j)\n® N \ij) = ipa z u ®S + <j+ ®W (28) 

where S = ( Ift ^- 1 ! V ) and W = i(4>\{S+ + 
S_)Q® N - X \ ip). Moreover, 

SS^ = \p\H(g)SS^-ipa-<S)SW^-(ip)*a+<Z>WS^+cr+a~(g)WWl 

(29) 

Now, on the one hand we see that the action of the left dissi- 
pator T> L leads to 

V L (S&) = 2T\p\ 2 a z u ®SS ] +Tip<j-®SW^ ' +T(ip)*<j+®W& . 

(30) 

On the other hand the left contribution of the unitary part of 
the Lindblad equation leads to 

i[H,S&]\\ L eft = -{ip+{i P Y)cj z u ®S^ 

-a- ® (l-2if L (ip))SW^ 
-<r+ <g> (l + 2if L (ip)*) WSHll) 

Comparing the two contributions gives the solution for the 
representation parameter 
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r - 2ip 



(32) 



The right boundary is treated along the same lines. The 
right-end state is polarized in the (x, z) plane in a direction v 
generated by the rotation 



U = e l 2 a 



(33) 



where we take as reference the — z-direction. Going through 
similar steps as above we impose the cancellation 



SI(0l) - i (S Z (0 L )< + S+{0 L )a+ + S-{B L )aZ) (24) 



S + (9 R M)=0. 



(34) 
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With (O this yields to 

i) = -tfm(6 R /2). (35) 

In order to fulfill the stationarity condition dT8b . together with 
( 1231 one needs to impose f R = —f L such that the repre- 
sentation parameter takes the value given in (f32l >. Interest- 
ingly, this condition turns out to allow for the inclusion of 
a Dzyaloshinsky-Moriya interaction in the XXZ-Hamiltonian 
Hill which is the key ingredient in the Lagrange-multiplier 
approach of I2H1 to current-carrying states of quantum spin 
systems. 

In conclusion, the solution of the completely polarized 
twisted case with a polarization on the left in the (y, z) plane 
and in the right in the [x, z) plane is given by the Matrix Prod- 
uct Ansatz for S in the form © with coherent state parameters 
d27l i. d35l l and representation parameter d32li. At 8l = Or = 
and vanishing boundary fields f R = f = one recovers the 
untwisted solution lfl2tl . 

The model with a twist is fundamentally different from the 
untwisted one, which can be seen by studying one- and two- 
point functions in the steady state. Note that in the isotropic 
model, all three spin projections <r x , a\ and a z n are locally 
conserved, i.e., j- t <r% = j«_ l n ~ j^ n+1 , where ft n+1 = 

2 Y J /3^ £ ap 1 (T ri (jl n+i i e <*Pi being Levi-Civita symbol). This 
leads to three different steady state currents j a = (]" n+1 ) 
for a = x,y,z. In the untwisted model (Or, = Or = 0) 
two out of three one-point correlations vanish in the steady 
state, (a x ) = ((7%) = for all n, corresponding to trivial 
flat x— and y— magnetization profiles along the chain. Also, 
two out of three spin currents are completely suppressed in 
the untwisted setup j x = j y = 0. In a model with a twist, 
neither of the one-point functions vanishes, and all three spin 
currents j x ,j y ,j z are generically nonzero. 

In order to see this, we note that in the usual untwisted 
model with Lindblad operators being creation/annihilation op- 
erators lfTll,[l2tl the steady state is invariant under a parity sym- 



metry p = UpU where U = (er z ) 0JV = U^ 1 [H. Any phys- 
ical observable that changes sign under the parity operation 
has to vanish in the steady state, e.g. (er*) = Tr (Ua x Up) = 
— Tr (a x p) = — (<r x ), from which (a*) = follows. In this 
way one readily obtains (af % ) = (cr^) = j x = j y = 0. 

In the isotropic model with a twist the parity symmetry is 
broken, but its place is taken by another symmetry, which we 
specify here for twisting angles 8^ = —6l = tt/2: It involves 
left-rightreflectioni?(A(g)5(g)...<g)C) = (C®....®B®A)R, 
global rotation in the XF-plane U ro t = diag(l,i)® N and 
v = an d reads p = VpVK where V = T, x U rot R 

Il20ll . It is straightforward to check that neither of the set of 
observables (a"),j a , changes sign under the V symmetry, 
and therefore they are generically nonzero. The symmetry 
V does, however, give rise to nontrivial relations between the 
observables, e.g. j x — —j y , (cr^) = — (a^_ n ), etc.. 

The major novelty of our approach is the exact MPA solu- 
tion by a quadratic algebra which turns out to be the symmetry 
algebra of the unitary evolution of the bulk part of the Hamil- 
tonian. Remarkably, this MPA solves the stationary Lind- 
blad equation even though both the quantum boundary fields 
and the Lindblad dissipators destroy this symmetry. We ex- 
pect that Lindblad equations for other open boundary-driven 
many-body quantum systems with a q-deformed non-Abelian 
bulk symmetry can be solved in a similar fashion. Since repre- 
sentations of the corresponding quantum algebras are known, 
exact results for observables become available. An open prob- 
lem is the relationship of the MPA to the integrability of the 
bulk Hamiltonian and hence to the extension of the MPA ap- 
proach to dynamical observables. Work on eigenfunctions 
ll22i l23ll and recent exact results by Eisler for the density ma- 
trix with bulk Lindblad terms j24ll hint at this possibility. 
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